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ABSTRACT 

Supernova remnants (SNRs) are believed to be the major contributors to Galactic cosmic rays. The 
detection of non-thermal emission from SNRs demonstrates the presence of energetic particles, but 
direct signatures of protons and other ions remain elusive. If these particles receive a sizeable fraction 
of the explosion energy, the morphological and spectral evolution of the SNR must be modified. To 
assess this, we run 3D hydrodynamic simulations of a remnant coupled with a non-linear acceleration 
model. We obtain the time-dependent evolution of the shocked structure, impacted by the Rayleigh- 
Taylor hydrodynamic instabilities at the contact discontinuity and by the back-reaction of particles at 
the forward shock. We then compute the progressive temperature equilibration and non-equilibrium 
ionization state of the plasma, and its thermal emission in each cell. This allows us to produce 
the first realistic synthetic maps of the projected X-ray emission from the SNR. Plasma conditions 
(temperature, ionization age) can vary widely over the projected surface of the SNR, especially between 
the ejecta and the ambient medium owing to their different composition. This demonstrates the need 
for spatially-resolved spectroscopy. We find that the integrated emission is reduced with particle back- 
reaction, with the effect being more significant for the highest photon energies. Therefore different 
energy bands, corresponding to different emitting elements, probe different levels of the impact of 
particle acceleration. Our work provides a framework for the interpretation of SNR observations with 
current X-ray missions (Chandra, XMM-Newton, Suzaku) and with upcoming X-ray missions (such 
as Astro-H). 

Subject headings: ISM: supernova remnants - Instabilities - ISM: cosmic rays - Acceleration of particles 
- Methods: numerical 



1. INTRODUCTION 



Young supernova remnants (SNRs) are a key link be- 
tween stars and the Galactic interstellar medium (ISM). 
Understanding their morphological and spectral evolu- 
tion is important to assess both explosion mod els of su- 
pernovae, o f which SNRs are fingerprints (e.g. |Badenes] 



et al.|2008 ) and large-scale mo dels of the turbulent IS 
for which SNRs are seeds (e.g. Joung fe Mac Low|2006[ ) 



A SNR is produced by the interaction of the material 
ejected during a supernova with the ambient medium 
(Chevalier 1977). The ejecta, being highly supersonic 
with typical Match numbers initially greater than 1000, 
induce a shock wave which propagates ahead of them in 
the surrounding medium. This shock, called the forward 
shock, accelerates and heats the ISM. As the ejecta push 
the layer of swept-up ISM, they get abruptly decelerated, 
inducing another shock wave that propagates backwards 
towards the center of the SNR. This shock, called the re- 
verse shock, decelerates and heats the ejecta. Although 
the basic dynamics of the shocked structure can be stud- 
ied with ID models, only 3D simulations can properly 
assess various important aspects such as asymmetric ex- 
plosions, mixing inside the ejecta, or interaction with a 
non-uniform ambient medium (in density or in magnetic 
field, for instance). In particular, the interface between 
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the shocked ejecta and the shocked ISM, called the con- 
tact discontinuity, is known to be hydrodynamically un- 
stable: because the dense ejecta are being decelerated 
by a medium of lower density, the interface is subject 
to the Rayleigh- Taylor instability, which prod uces a dis 
tinctiye pattern of fingers and holes (see e.g. |Chevalier 



er 
1 



eTliin|T992l IDwark adas 2000J |Wang fe Chevalier||^00' 
Blondm & Eliison||200ip. It is important to know how 
such structures will affect images of these objects, which 
are only seen in projection, and how this may affect their 
diagnostics. 

In addition to their key role in the dispersion of heavy 
elements synthesized inside stars and during supernovae, 
SNRs are widely believed to be the main production sites 
of Galactic cosmic-rays up to energies of about 10 15 eV. 
A fraction of the particles swept-up by the blast wave 
is expected to gain energy through the diffusive shock 
acceleration mechanism (DSA), a first-order Fermi pro- 
cess which relies on both the velocity discontinuity at the 
shock and the presence of magnetic turbulence on both 
sides of the shock. Whether the same mechanism could 
also operate at the re yerse shock propaga ting back into 
the ejecta is debated (Ellison et al. 2005). In any case, 
growing evidence has accumulated in the last decade 
ma king the case fo r SN Rs be i ng pa rticle accelerators; 
see 
on 



Reynolds (20081) and Vink (20121 for rec ent reviews 
he high- ener gy emission from SNRs, and|Ferrand fc| 



Safi-Harb ( |2012[ ) for an up-to-date catalogue of X-ray 
and 7-ray observations of Galactic SNRs. However, a 
direct proof of the acceleration of hadrons in these ob- 
jects is still lacking. Particles detected without ambi- 
guity are highly energetic electrons, which radiate far 
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more efficiently than protons. In the radio regime, the 
synchrotron radiation from electrons has been routinely 
detected for half a century, and is in fact a cornerstone 
for the identification of SNRs. In X-rays, the tail of this 
same radiation has been detected in several objects in 
the last 20 years, and was even found to dominate the 
emission in a few SNRs. In gamma-rays, several rem- 
nants have finally been detected in the last decade; the 
pending issue remains to disentangle the leptonic and 
hadronic contributions associated with I nverse Comp ton 
scattering and pion decay, respectively ( |Gabici|2008 ). It 
is thus important to be able to simulate trie non-thermal 
emission from SNRs . Su ch work has been und ertaken by 
Petruk et al. (2009b I and Orlando et al. (2011 1, assuming 
a given distribution of accelerated electrons on top of 3D 
MHD simulations of the remnant from SN1006. Using 



a model of acceleration, Lee et al. ( 20081) computed the 
broad-band emission from both protons and electrons, on 
top of a ID hydrodynamic model of the remnant's expan- 
sion. At this point it is worth noting that, if SNRs accel- 
erate particles as efficiently as we think, then this must 
affect their own dynamics in a sizeable way. If a sub- 
stantial fraction of the explosion energy is diverted into 
the acceleration of particles at the shock, the evolution 
of the latter must substantia l ly deviate from usual hy - 
drodynamics (Chevalier 1983 Decourchelle et al. 2000). 
This, in turn, should have a visible impact on the ther- 
mal emission from the remnant. Therefore, probing the 
morphological and spectral evolution of the SNR itself 
provides an indirect, but promising, way to assess the ef- 
ficiency of particle acceleration in young remnants and, 
importantly, the acceleration of hadrons given that lep- 
tons have a negligible inertia. Indeed observational ev- 
idence for modified mor phologies has been r eported in 
the remnants of Tycho (|Warren et~aLl|20 05), SN 1006 
( ICassam-Chenai et al.|2008[|Miceh et al.| 1 2009f and Cas A 
( Patnaude &: Fescn 2009J . Thus a few teams have begun 
to address the effect of back-reaction of particles on the 
evolu tion of a SNR and i ts impact on the thermal emis- 
sion. |Ellison et al. (2007) investigated both the thermal 
and non-thermal emission from a spherically sy mmet- 
ric SNR, using a semi-analytical model for DSA ( |Blasi| 
et al. 2005[ ) and an approximate treatment for the ion- 
ization state of the plasma. Patnaude et al (2009 20101 



expanded on this work by making a direct coupling of 
an ionization and emission code with a ID hydrodynam- 
ical code, to study how acceleration impacts the shocked 
ISM behind the forward shock (the emission fr om the 
eject a was not considered). Using this model, Castro 



et al. (2011) produced broad-band spectra with both" the 



thermal and non-thermal emission. iFerrand et al.|(|2010|) 
coupled the same kind of acceleration model to a 3D 
code, enabling the study of the development of the in- 
stabilities over the edge of the ejecta concomitantly with 
the shrinking of the shocked region caused by energetic 
particles. The hydrodynamic results show that the posi- 
tions of the shock waves in Tycho's SNR can readily be 
explained by the presence of energetic ions. 

Here we extend this work, by computing the emis- 
sion from SNRs undergoing efficient particle acceleration. 
In the present paper we compute the thermal emission 
(the non-thermal emission will be presented in a sub- 
sequent paper). More precisely, we concentrate on the 
X-ray emission, which is quantitatively the largest emis- 



sion coming from the hot interior of SNRs, and which 
is qualitatively an important source of diagnostics for 
these objects. Thermal X-ray emission allows for the 
derivation of global parameters of the remnant, such as 
the density and expansion velocity. Comparison of the 
observations with models and simulations also allows the 
assessment of key parameters of the acceleration process, 
such as injection fraction or energetics efficiency. 

The outline of the paper is as follows. In Section [2j 
we present our method. We recall how we perform 3D 
hydrodynamic simulations of the remnant evolution, cou- 
pled with a model of non-linear acceleration. We explain 
how these simulations are post-processed to obtain the 
non-equilibrium ionization state of the plasma and com- 
pute its emission. In Section [31 we present our results. 
We first show some useful hydrodynamic and thermo- 
dynamic quantities, then discuss synthetic spectra and 
projected maps from our model SNR. In Section [4j we 
summarize our main results in link with observations, 
and discuss the main physical assumptions underlying 
these. In Section [5] we present our conclusions and per- 
spectives. 

2. MODEL 

In this section we present our 3D numerical model, 
which computes the mor pholo gical and spectral evolution 
of a young SNR (Section 2.1 1, including the acceleration 
of particles at the forward shock and the back-reaction 
effects (Section |2~2j). 



2.1. Morphological and spectral evolution of a young 

SNR 

In our approach computations are done in two steps: 
first we comp ute the hydrodynamic evolution of the rem- 
nant (Section 2.1.1 1 with some assumptions on the ther- 
modynamics of the shocked material (Section |2.1.2 ), then 
we compu te the ionization state of the pl asma (Sec- 
tion 2.1.3) and its thermal emission (Section 2.1.4[ ). 



2.1.1. Hydrodynamics 

We consider that the SNR is initially spherically sym- 
metric. The structure of the shocked region then depends 
on the radial density profiles of the ejecta (p e jc) and of 
the ISM (pism)- We assume power-law profiles in both 
media: 

p cjo (r,t)ocr-"^ 3 , (1) 



Pism{r,t) oc r~ 



(2) 



so that the radial profiles of all hydrodynamic quantities 
(density, velocity, pressure) are initially self-similar, and 
can be computed by the integra tion of ordinary differ- 
ential equations (Chevalier 1982 1. Given equations ([I]) 
and ([2]), dimensional analysis shows that the radius of 
the contact discontinuity increases with time as 



a(t) =i(»-3)/(«-s) ) 



(3) 



These solutions are valid until the reverse shock reaches 
the central core of the ejecta. After the initial ejecta- 
dominated phase, the SNR enters another phase, still 
adiabatic but now dominated by the dynamics of the 
swept-up ISM, which is commonly called the Sedov 
phase. To track the evolution of the SNR over time, with- 
out any limiting assumptions, we make time-dependent 
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numerical simulations, where the profiles are initialized 
with the self-similar solutions and then evolve freely. We 
use the code RAMSES, which contains a Eulerian hy- 
drodyna mical solver ba sed on a second-order Godunov 
method ( |Teyssier||2052| ). We run it with the MUSCL 
scheme and an exact Riemann solver. The physical do- 
main is discretized on a Cartesian grid. As we want to 
catch the development of the Ray leigh- Taylor instabili- 
ties, we use a 3D grid. As the SNR can expand by sev- 
eral orders of magnitudes, and most interesting features 
are present inside the shocked region, we use a comoving 
grid, in which the system of waves is quasi-stationary. To 
obtain this, we perform a change of variables based on 
appropriate powers of the scale factor defined by equa- 
tion pi), and we add a non-ine rtial force to the Eule r 
equations (for more details see Fraschetti et al. 2010). 
The resolution of the grid around the shocks is further 
increased by using RAMSES tree-based adaptive mesh 
refinement (AMR). 

2.1.2. Thermodynamics 

The shocked plasma (both the ISM downstream of the 
forward shock, and the ejecta downstream of the reverse 
shock) is heated to such high temperatures that it shines 
in X-rays. In the reference case without any losses, that 
is, as long as a negligible amount of energy is channeled 
into energetic particles or lost through radiationFjthe av- 
erage temperature immediately downstream of the shock 
(Ts) depends solely on the shock speed (ug) according to: 



we use the m ethod of Co x fc A nderson (1982), based on 



rr 3 m 

Iq = V 

b 16fc B 



2 x 10 7 K 



1000 km/s 



(4) 



where m = jj,m p is the mean molecular weight (m p being 
the proton mass and fj, the weight factor accounting for 
the composition of the plasma) and ks is the Boltzmann 
constant. 

The emissivity of the plasma depends on the temper- 
ature of the electrons (T e ) which excite and ionize ions. 
It is useful to introduce the ratio /3 of the electrons and 
protons temperatures: 



(3 



T 
1 p 



(5) 



which is a function of time and position. The value of 
this ratio at the shock front (/3s) is poorly unknown, 
and is therefore a free parameter of the model. Theo- 
retically, it could be anywhere between the masses ratio 
(/3s = m e /m p ) and full equilibration (/3s = 1). Ob- 
servations suggest prompt equilibration for shock speeds 
below about 400 km/s, but intermediate degrees of equi- 
libration for faster shocks, with /3s ten ding towards zero 



rough ly as v s (see Rakowski 



2005 



and 



Ghavamian et al. 



2007 for reviews). Downstream of the shock, electrons 



will be heated through Coulombian interactions with 
ions, so that T e will progressively equilibrate with T p j^] 
To compute T e in any given cell in the shocked region, 

5 Note that although in this paper we are computing the ra- 
diation from the remnant, we restrict our study to the so-called 
non-radiative phases of its evolution, meaning that radiation has 
no impact yet on its dynamics. For cosmic-rays the situation is re- 
versed: they are thought to be mostly produced during these early 
phases, and thus are fully accounted for in our approach. 

6 Assuming that ions species are heated according to their mass, 



the model of |H pTcTT 1979 . This moW assumes an 
adiabatic evolution ot the plasma. It provides the cur- 
rent value of the temperature ratio /3, given its value at 
the shock /3s, the current density n and temperature T 
of the cell, and its age t§ defined as 



7fc(x,t) = f dt' 



(6) 



where ts is the time when the material considered was 
shocked. In the numerical simulations, we keep track of 
rs with the inclusion of a passive scalar in RAMSES, 
which is incremented in the shocked region only. 

2.1.3. Non-equilibrium ionization 

Free electrons in the plasma excite and ionize the 
atoms, until an equilibrium is reached between ioniza- 
tion and recombination for all spec ies. We use ioniza- 
tion and r e combination rates from [Arnaud fc Rothen 

and [Mazzotta 



, |Arna 

et al. ( |l998[ r The final ionization state depends only 



Aug (1985)JArnaud fc Raymond ( 1992 1 



on the electron temperature (T e ), but remnants we con- 
sider here are much too young to have reached ioniza- 
tion equilibrium, so that the ionization state also de- 
pends on the electron density (n e ) and on the time (i) 
elapsed since the plasma has been shocked. Ionization 
timescales for var ious elements, as a func tion of temper- 



ature, are given in Smith & Hughes (2010). It is apparent 



that no species heavier than Carbon can be ionized for 
n e x t < 10 7 cm~ 3 .s, and that equilibrium of all ions is 
granted for n e x t > 10 13 cm~ 3 .s. In between, the ac- 
tual ionization state must be computed. The value of 
the electron density, n e , at any time, t, depends on the 
ionization state of the plasma, so that the computation 
of its ionization state should in principle be coupled with 
its hydrodynamic evolution at each time step. We sim- 
plify the problem by deferring the computation of the 
actual n e to the end of the simulation, recording only 
the evolution of the total density over time. The hydro- 
dynamic code natively works with the mass density, p, so 
we added another passive scalar, tt, in RAMSES defined 
as: 



n(x,t)= f P (x,t').dt' 



(7) 



where we recall that ts is the time when the material 
considered was shocked (note that we use the convention 
of denoting quantities defined since the shock crossing 
with Greek letters). So after material is shocked, we 
keep track of its ionization age as it is advected down- 
stream, and at the end of the simulation we compute the 
non-equilibrium ionization state in each cell. To do so 
we use the so-call ed exponentia tion method, introduced 
indepe ndently by Masai ( 1 984 ) and | Hugh es fc Helfand| 
(1985), and reviewed in Smith & Hughes (2010 ). In gen- 
eral, the ionization evolution ot an element of atomic 
number Z is given by a coupled set of Z + 1 equations; 

and given that their various masses are not widely different (com- 
pared to the electron mass), we consider that all ion species quickly 
get roughly the same temperature after crossing the shock front, 
and focus only on the progressive equilibration of electrons with 
ions. 
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in this method, these are recast as Z uncoupled equa- 
tions, after diagonalization of the rates matrix. These 
Z first-order differential equations are integrated to ob- 
tain the ionization fractions of the element, for a given 
electronic temperature (T c ) and ionization age (J n c dt). 
The ionization age is defined implicitly, as n c depends 
on the ionization fractions, which can only be computed 
for some value of J n e dt, hence the need for an itera- 
tive computation. A first guess for the ionization age 
is obtained from the tracer t\, assuming that all species 
are ionized once - an assumption that is sufficient for the 
shocked ISM, where Hydrogen dominates, but not for the 
shocked ejecta, where heavy elements are not negligible. 
This value is used to compute the ionization fractions, 
from which the actual value of n e is obtained. Then the 
ionization fractions are re-computed, and the thermal 
emission is computed. 

2.1.4. X-ray thermal emission 

The hot ionized plasma located between the two shocks 
of a young SNR is a source of photons with energies of the 
order of keV, which is now routinely o bserved with space- 
based X-ray instruments (see e.g. 



Decourchelle 



2005 



for a review). This X-ray emission is a function of the 
electronic density and temperature of the plasma, and 
of the ionization fractions of each species (which them- 
selves depend on density and temperature as explained 
in the previous section). We compute this thermal emis- 
sion as a function of the photon e nergy with a spectral 
code from Mewe et al. (1985 1986). The code computes 
both the continuum emission (tree-free, free-bound and 
two-photons radiations) and the emission lines (from ex- 
citation and recombination processes), for each of the 
following 15 elements: Hydrogen (H, Z = 1), Helium 
(He, Z = 2), Carbon (C, Z = 6), Nitrogen (N, Z = 7), 
Oxygen (O, Z = 8), Neon (Ne, Z = 10), Sodium (Na, 
Z = 11), Magnesium (Mg, Z = 12), Aluminium (Al, 
Z = 13), Silicon (Si, Z = 14), Sulfur (S, Z = 16), Argon 
(Ar, Z = 18), Calcium (Ca, Z = 20), Iron (Fe, Z = 26), 
and Nickel (Ni, Z = 28). 

The total thermal emission from any given cell is the 
sum of the emission from the ejecta and from the ISM. 
The distribution of the ejecta is known thanks to another 
scalar variable passively advected with the fluid density, 
the fraction of ejecta /, defined so that the ejecta density 
is / x p and the ISM density is (1 — /) x p. The thermal 
emission of species i is proportional to n e x n„ with n e 
the electrons density and rii the ion density, which are 
computed from p according to the abundances and ionic 
fractions. The ejecta and the ISM are distinguished by 
their composition: we assume solar abundances for the 
ISM (we use standard values from Anders & Grevesse 
1989 1, and metal-en hanced abundances tor the ejecta (we 

For reference, 



2010 



adopt the model of Maeda et al. 

this corresponds to molecular weights oT/i ~ 1.9 in the 
ejecta and fi ~ 0.6 in the ISM, if we assumed full ioniza- 
tion (the most abundant, light species in the ISM are in- 
deed likely to be fully ionized; the most abundant, heavy 
species in the ejecta are not likely to be). Note that, 

7 Other composition models are available for the ejecta (see for 
instance the reference work bv|Nomoto et al.|1984| and a more re- 
cent work by Ba denes et al.|2005b[ ), but comparing them is beyond 
the scope of this paper. 



in each medium, we do not consider spatial variations 
in the abundances, as in this first study of the thermal 
emission from 3D simulations, we are mostly interested 
in disentangling projection effects, independently of the 
supernova which gave birth to the structure. 

The emission is computed in each shell of the shocked 
region between the reverse and the forward shocks. It is 
then projected along one axis of the cartesian grid, prop- 
erly taking into account the tree-based AMR structure. 
We thus obtain 2D synthetic maps of the emission of the 
remnant. Note that, even though our final product is 2D, 
it is mandatory to make 3D simulations to properly cap- 
ture the development of the Rayleigh- Taylor instabilities 
at the contact discontinuity. This is all the more im- 
portant since the thermal emission is dominated by the 
shocked ejecta which exhibit a complicated morphology. 

2.2. Acceleration and back-reaction of particles 

Strong and magnetized shocks, such as the ones en- 
countered in young SNRs, can accelerate particle to high 
energies. It is not the purpose of this paper to study the 
fate of these particles in detail, and in particular their 
non-thermal emission which will be presented in a sub- 
sequent paper. Our point here is to show how the thermal 
emission from the SNR is affected by particle accelera- 
tion. In this section we exp lain h ow we model acceler- 
ation at the shock (Section 2.2.1), how we handle the 
bac k-react ion of accelerated particles on the flow (Sec- 
tion 2.2.2), and how this affects the diagn ostics needed 
to compute the thermal emission (Section 2.2.3). Note 
that, although our model is general, we choose to use it 
at the forward shock only, because it is unclear whether 
the reverse shock can be also a n effi cient accelerator (see 
more in the discussion, Section 4.5). 



2.2.1. Diffusive shock acceleration 

A fraction of the particles swept-up by the blast wave 
is expected to be accelera ted through t he diffusive shock 
accele ration mechanism (Drury 1983 Malkov & Drury 
20011. In the test-particle regime, this mechanism pro- 
duces a power-law energy spectrum with an index (s) 
that is universal, because it is entirely determined by 
the compression ratio (r) of the shock, which for strong 
shocks in the ISM is always r — 4. However, as the pres- 
sure of particles builds up upstream of the shock, the gas 
is pre-accelerated, hence the formation of a precursor. 
As the distance a particle can travel ahead of the shock 
is a function of its energy, particles of different energies 
sample different parts of the precursor, where they expe- 
rience different velocity jumps, and therefore the slope s 
becomes energy-dependent: spectra get concave. Solving 
explicitly the coupled system of the fluid and the parti- 
cles as a function of time can be done numeric ally (e.g. 
Ferrand et al.|2008||Zi"rakashvili fe Ptuskin|20"TT| ), but the 
high numerical cost limits the possibility to run 3D mod- 
els of realistic SNRs. Instead, w e use t he sem i-analytical 
model of acceleration of Blasi ( |2002| |2004[ ). Given the 
shock properties (velocity and upstream conditions) , this 
model jointly solves the particle spectrum and the fluid 
velocity profile as functions of the momentum of parti- 
cles. Note that it does not solve the time evolution of 
the system, but only finds quasi-stationary solutions, so 
that we have to rerun it after each hydro time-step. This 
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procedure is justified at most energies, but breaks down 
close to the highest energies when the acceleration time 
gets of the same order as the age of the SNR (the ac- 
celeration time is comp uted with the classical formula 
derived by Drury ( 1983[ ) in the linear regime, but using 
the modified velocities upstream and downstream of the 
shock to account for non- linear effects). The maximum 
energy can also be limited by the lack of geometrical con- 
finement: the model includes self-consistently the escape 
of particles upstream o f the shock, which carry energy 
away from the system (Blasi et al. 2005). Assuming a 
uniform explosion in a uniform ISM, we expect the shock 
properties to be symmetric, and so at any given time we 
run the acceleration model only once, using azimuthally 
averaged shock properties. The shock fronts are identi- 
fied in the hydro simulation through their strong pressure 
gradients. 

Even though the magnetic field does not play a dy- 
namical role in the evolution of the SNR (and is thus not 
activated inside RAMSES), the nature and the level of 
magnetic turbulence in the vicinity of the shock is cru- 
cial to particle acceleration. Both the age-limited and 
size-limited maximum energies of particles are computed 
for a given diffusion coefficient D. We make the common 
assumption of Bohm diffusion, that is 

D(p) = Do-? with D = 3 - 1 ° 2 : C n' S ' 1 PS 
y/l +p 2 S 0»G) 

where p is the momentum of the particles in m p c units 
and B is the intensity of the magnetic field. Both theo- 
retical and observational evidence suggest that energetic 
particles c an substantially am plify B upstream of the 
shock (see Bykov et al. 2012 for a review), but we do 
not include this effect in our model at this point because 
of the large uncertainties and complexities in the mecha- 
nisms at work; we will address this issue in a subsequent 
paper where we comput e the non- t herma l emission of 
accelerated particles (see Lee et al. (2012) for other re- 
cent numerical developments on this issue). However, we 



here include t urbulent heating in the precursor (Amato 
& Blasi 2006): we assume that Alfven waves are effi- 
ciently damped in the precursor, so that the gas gets 
compressed non-adiabatically. This tempers the inten- 
sity of the back-reaction, so that the particle pressure 
does not get unreasonably highj^] 

2.2.2. Particle back-reaction 

When it comes to the SNR evolution, what matters 
foremost is the modification of the shock compression. 
As a precursor is formed]^] the shock discontinuity is re- 

8 If, instead, we had assumed that amplified Alfven waves are 
advected to the shock, the jump relations would be modi fied, pro- 
viding ano ther kind of regulation of the shock structure ||Caprioli| 
|et al.|2008| |. The de gree of modification of the SNR would be ot the 
same order, although a little higher. As far as the thermal emission 
is concerned, choosing a different recipe for magnetic field amplifi- 
cation would amount to choosing a slightly different value for the 
injection fraction (for the effec t of n on the shock modification, see 
Fig. 2 oFFerrand et al. 2010). Implications for the non-thermal 
emission will be discussed m a subsequent paper. 

9 In our hydrodynamic simulations we do not explicitly handle 
this precursor upstream of the shock. Its impact on the X-ray 
emission is likely very small. It can however be detected in the 
optical domain as if affects the Balmer lines of Hydrogen and the 
forbidden lines of heavier elements (e.g. |Rakowski et al.|20l"T} . 



duced to what is called the sub-shock, of compression 
ratio r su b < 4, whereas the overall compression (from far 
upstream to downstream) is increased, with a compres- 
sion ratio r to t > 4. To take that effect into account in 
our simulations, we use a pseudo-flui d with variable adi 



abatic index 7, as we already did in |Ferrand et al. (20101 



(we recall that 7 is defined so that 

P 



7-1' 



(9) 
(10) 



where p, P, e are the fluid's mass density, pressure and 
internal energy, respectively). In our version of RAM- 
SES, 7 is a variable of each cell, and is modified in order 
to account for the presence of energetic particles. The 
pseudo-fluid thus represents both the thermal (shocked 
fluid) and non-thermal (accelerated particles) popula- 
tions. To enforce the back-reaction of accelerated par- 
ticles on the SNR structure, we thus proceed as follows: 
at each time-step, we compute the index (jcs) which will 
produce the same ratio (i?tot) as predicted by the accel- 
eration model, and affect it to the cells located just up- 
stream of the shock front. 7 e g thus drops down from the 
fiducial value of 5/3, close to the value of 4/3 appropri- 
ate for relativistic particles, or even slightly below that 
due to the escape of particlespj Then j e g is advected 
downstream of the shock, inside the shocked region, so 
that each fluid element remembers modifications induced 
by particle acceleration at the time it was shocked. This 
modification of 7 makes the fluid more compressibl e, so 
that the shocked regi on gets denser and narrower (De- 
courchelle et al 



2000) 



The importance of back- reaction effects obviously de- 
pends on the amount of energy channeled into energetic 
particles, or alternatively on the fraction of particles that 
enter the acceleration process at the shock front. The lat- 
ter mechanism, known as the injection, is unfortunately 
poorly known, and therefore often paramet rizedF^] We 
adopt here the thermal leakage recipe of ( |Blasi et al. 
2005 1 and set the fraction rj of particles crossing the shock 
entering the acceleration process to be 



V 



Sir 1 / 2 



(r- su b-lK 3 exp(-£ 2 ) 



(11) 



where r su b is the compression ratio at the (sub-)shock 
(see next section) and £ is a free parameter parametrizing 
the energy of swept-up particles, defined as 



Pinj 

Pth 



(12) 



where p ln j is the injection momentum of particles and 
Pth is the mean thermal momentum immediately down- 
stream of the shock. £ is expected to be in the range 

10 A more rigorous treatment would be to split the two effects: 
lower 7 e ff to account for the presence of relativistic particles, add 
an energy sink to account for the escape of some of these particles. 
This was done by |Patnaude et al.] ( |2009} , who observed that such 
changes are important lor tne later stages of the SNR evolution, 
but not at the age considered here. A significant improvement 
over our current model would be a multi-fluid approach, where 
each population (thermal and non-thermal) has its own 7. 

11 It is however not the case in Monte-Carlo simulations of the 
kind of |Ellison fc Eichler| p984) . 
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2-4. The fraction r\ obtained from Eq. [TTJ is applied to 
the total mass flux crossing the shock, regardless of its 
detailed composition: the acceleration model only asks 
for the proton number density, which is a reasonable sim- 
plification for the ISM swept-up by the forward shockFj 
Finally, note that back-reaction effects might be al^ 
ready significant at the small age at which we initialize 
the simulation using self-similar profiles. The pressure of 
rclativistic particles, computed from the same accelera- 
tion model, is therefore includ ed in the computation of 
these profiles (Chevalier 1983). 



2.2.3. Properties of the modified plasma 

As already explained, the emission at each point de- 
pends on the electrons density (n e ) and temperature 
(T e ). These quantities are derived from the mass den- 
sity (p) and the protons temperature (T p ) of the hydro 
code as explained before. But because of our use of a 
pseudo-fluid to implement the back-reaction of particles, 
we need to extract the properties of the thermal com- 
ponent only. The density of the pseudo-fluid is within a 
very good approximation the same as the density of the 
thermal fluid, given the low fraction of particles acceler- 
ated at the shock. Regarding temperature, a distinction 
must be made between the ejecta and the ISM. In the 
shocked ejecta, the pressure of the pseudo-fluid is actu- 
ally the pressure of the real fluid, given that we do not 
consider acceleration at the reverse shock. Hence the 
temperature in the shocked ejecta (T e j) can simply be 
obtained as: 



/im p P 
p k B 



(13) 



In the shocked ISM, the pressure of the pseudo-fluid con- 
tains an important contribution from energetic particles 
- in fact, it can easily be dominated by it - so that we 
cannot use the same method. We determine the tem- 
perature in the shocked ISM (Tism) a posteriori, using 
adiabatic evolution downstream of the shock: 



- ISM 



T s x 



7-1 



(14) 



where T = T (r, t) and n = n (r, t) are the temperature 
and density of a cell at radius r and time t, and Ts = 
T(rs,ts) and ng = n(rs,is) are the temperature and 
density immediately downstream of the forward shock 
at radius rg at the time it was shocked t$. We recall 
that each cell knows how much time rs elapsed since 
its material was shocked (Eq. [6]). We also need to keep 
track of the history of shock modifications: the values 
of the density ns and temperature Tg just downstream 
of the shock, obtained from the acceleration model, are 
recorded at each time step. 

Knowing the hydrodynamic and thermodynamic prop- 
erties of the thermal plasma, we can compute its ioniza- 
tion state and thermal emission as explained before, now 
taking into account the back-reaction of accelerated par- 
ticles on the SNR. 



regarding ' 
process. 



See |CaprioIi et al.|(|201l| for recent development 
ding the impact ol heavy species on the diffusiv 



2.3. Parameters 

Our goal in this paper is not to make a detailed mod- 
eling of a specific object, but rather to point out the im- 
portant physical aspects to consider when interpreting 
observations. So we do not to explore the full parameter 
space, but rather use a single set of represent ative pa- 
rameters (a s in our preliminary work shown in Ferrand 
eFaI1[2TTro| . 



s of the model 
diffusive acceleration 



We are considering the remnant of a type la event, 
and are using Tycho's SNR as a reference, but we are 
not attempting a fit of this object. We set the explo- 
sion energy to E = 10 51 erg, and the ejecta mass to 
M j = 1.4 solar masses, and adopt a power-law pro- 
file with index n = 7. We assume a uniform ambient 
medium (that is, a power-law profile with index s = 0) 
of rather low density n,H,o = 0.1 cm -3 . We use these pro- 
files to initialize the hydrodynamic code at a small age 
of t = 10 years. We then let the code evolve the profiles 
until t = 500 years. Particles are accelerated continu- 
ously in time at the forward shock, for the injection level 
we take the common choice £ = 3.5, which typically pro- 
duces injection fractions of the order of a few 10~ 4 for 
the remnant considered. The injection level is expected 
to vary with the obliquity angle between the mean mag- 
netic field and the shock normal, but whether parallel 
or perpendicular shocks are more efficient accelerators 
is still debated, and the geometry of the magnetic field 
around the blast wave is not well known anyway. So we 
do not consider here the obliquity dependence of the in- 
jection, and instead use this fraction as an average value 
for the whole shock surface, at any given time. This as- 
sumption of overall symmetry seems reasonable for an 
object like Tycho (for a clearly bipolar object like SN 
1006 we refer the reader to the work by |Orlando et al. 
2007| and jPetruk et al. ||2009c|a[ ) . As time increases, "tne 
injection fraction is lowered as r sub is reduced, see Eq.[TT} 
We assume an ambient magnetic field B — 5 /iG for tne 
diffusion, leading to maximum energies of a few 10 5 m p c 2 
for protons. 

Hydrodynamic results are processed at the final age 
to compute the thermodynamic and ionization state of 
the plasma in each cell. The only free parameter intro- 
duced at this stage is the ratio /3s of electron and proton 
temperatures at the shock front, which for the sake of 
simplicity we fix here to /3s = m e /nip ~ 5. 5 x 10~ 4 (in 
the case of Tycho, observations reported by |Ghavamian 
et al.|[2001 suggest /3s < 20% at the forward shock). 

We then compute the thermal emission in each cell. 
There is no parameter to tune in the emission code. We 
just have to define the energy grid of photons, which we 
choose to be logarithmic, with 100 points per decade over 
2 decades in energy, from 0.1 keV to 10 keV - which is 
wide enough to contain most of the X-ray emission ob- 
served by current instruments. We recall that the emis- 
sion is computed only in the shocked region, as bounded 
by the (average) positions of the reverse and forward 
shock, and then projected along the line of sight. Spa- 
tially, we use an AMR grid with a formal resolution of 
1024 3 , so that we produce emission maps of 1024 2 pixels 
in each energy bin. The effect of numerical resolution is 
shown in Appendix [AJ Note that we only simulate one 
octant of the remnant, and therefore obtain projected 
maps of one quarter of the remnant. 
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3. RESULTS 

In this section we present a variety of synthetic emis- 
sion maps and spectra generated from our 3D simula- 
tions, in order to investigate the effects of acceleration 
on observational diagnostics of the SNR. In all the plots, 
we compare data obtained without including the back- 
reaction of accelerated particles (cases labeled 'OFF'), 
that is when forcing the (test-particle) linear regime, so 
that the shock structure is un-modificd; and data ob- 
tained including the back-reaction of accelerated parti- 
cles (cases labeled 'ON'), that is when letting the non- 
linear regime develop, so that the shock structure is mod- 
ified. We first show th e qua ntities on which the thermal 
emission relies (Se ction 3.1 ), before showing the emission 
itself (Section 3.2 1. 



3.1. Plasma state 

Here we show the quantities, diagnosed from the RAM 
SES code, needed to co mpute the emission: hydrody 
namic quant ities (S ection 3.1.1 1, thermodynamic quanti 



ties (Section |3.1.2 ), and ionization state (Section 3.1.3) 
We present all these data in the form of slices (in a plane 
z = perpendicular to the line of sight) in order to show 
the inner structure of the remnant. 

3.1.1. Hydrodynamics 

The mass density p is shown at the top panel of Fig- 
ure [II It ranges over one order of magnitude, from 
roughly 0.1 to 1 m^cm -3 . The fluid consists of different 
species, the density shown here is the sum of all nucleons 
over all elements. The fluid also consists of two media, 
the ejecta in the inner region and the ISM in the outer 
region. In between the two we observe the development 
of the Ray leigh- Taylor instabilities: fingers of ejecta are 
protruding inside the ISM, with mushroom-shaped heads 
that sometimes get disrupted. The contact discontinu- 
ity between the ejecta and the ISM gets a very irregular 
shape, and is often ill-defined as the two media get mixed. 
The bottom panel of the same Figure [l] sho ws the ioniza- 
tion age / n e .dt of the material (see Section 2.1.3 ). In the 
ISM, this age mostly reflects the time since the material 
was shocked, rising progressively from the forward shock 
to the contact discontinuity, which is highlighted by a 
sudden drop in density. The interface is much distorted 
though, so that old material can be found in a wide range 
of radii within the shocked region - a feature which makes 
3D simulations mandatory. In the ejecta, the ionization 
age remains lower, even though they are denser, because 
heavier elements take more time to get ionized after they 
get shocked by the reverse shock. With our parameters, 
the maximum ionization age is ~ 2.5 x 10 9 cm _3 .s in the 
un- modifie d case. According to the discussion of Sec- 
tion |2.1.3[ it is apparent that the plasma must be in a 
non-equilibrium state. 

The effect of acceleration can be seen when comparing 
the left and right sides of the figure: as part of the energy 
of the system is diverted into the acceleration of parti- 
cles, the shocked region shrinks. This means two things. 
First, it gets narrower in radius, here by about one t hird, 
which already m akes an observational diagnostics (Fer- 
rand et aT]|2010 ). Second, it gets denser close to the for- 
ward shock (where acceleration happens), here by a fac- 
tor of about 2. Whereas in the un-modified case the high- 



est densities are found in the ejecta, close to the reverse 
shock, in the modified case the same kind of values are 
obtained next to the reverse and forward shocks. The re- 
gion immediately downstream of the forward shock looks 
irregular when back-reaction is activated, but we warn 
the reader that at least part of it is caused by numerical 
instabilities triggered at the shock front. These insta- 
bilities, commonly known as carbuncle, develop because 
of th e quasi-stationarity of the shock in the comoving 
grid ( |Fraschetti et al.| 2010 ) , and are exacerbated when 
the shock is modified by particles^] The overall increase 
in density close to the shock, however, is perfectly real. 
And part of these post-shock turbulent features might in 
fact be real, and related to Rayleigh- Taylor structures, as 
the radial density profile is substantially modified in the 
shocked ISM. At any rate, the exact density profile close 
to the forward shock should not much affect the maps 
of thermal emission presented later, as long as the ejecta 
dominate, which is usually the case, and even more so 
when including back-reaction, so that we can pursue our 
study. As we do not include acceleration at the reverse 
shock, the region between the reverse shock and the con- 
tact discontinuity is not much affected by acceleration; 
however, it seems that the ejecta cannot protrude as far 
away in the ISM, as they approach the closer and denser 
outer shock edge. Because the shocked ISM is more com- 
pressed, the ionization age gets higher here, by almost a 
factor of 2. 

3.1.2. Thermodynamics 

Temperatures are show n on F igure [2] for protons and 
for electrons (see Section 2.1.2). As expected, the tem- 
perature is higher in the center of the shocked region, 
close to the contact discontinuity, where material was 
shocked at an earlier age (compare with bottom of Fig- 
ure [lj, when the shock speed was much higher. Over 
the course of the simulation, the speed of the blast wave 
in the ISM falls down from ~ 30, 000 km.s -1 at the ini- 
tial time t = 10 yr to ~ 5000 km.s -1 at the final time 
t = 500 yr. In the un-modified case, theprotons temper- 
ature is of the order of 10 9 K (see Eq. H. The electrons 
temperature (the one that matters regarding thermal X- 
ray emission) is much lower (by more than two orders of 
magnitude), given our assumptions that heating at the 
shock depends on the mass ratio and equilibration with 
protons is progressive. T c can still reach about 20 x 10 6 K, 
which is enough to produce intense X-ray emission. 

The effect of acceleration is important for the ISM, 
given our assumption of efficient acceleration at the for- 
ward shock. Because a significant part of the energy goes 
into acceleration of particles, the temperature of protons 
T p is much lower in the shocked ISM with back-reaction: 



by typically a factor of 4 in th e case simulated here 
effect is well known, see e.g. |Decourchelle et al 



This 



for an early ID modeling and Hughes et al. 



2000) 
20001) for a 



possible observational diagnostic. 'l p is also indirectly re- 
duced in the shocked ejecta, although by usually less than 

13 Although unphysical features could be smoothed with vari- 
ous numerical strategies (including using more diffusive Riemann 
solvers, or using a wobbling comoving grid), it appears difficult to 
completely avoid them while maintaining a sufficient resolution of 
the physical structures inside the shocked region. As a result, our 
current numerical setup does not allow us to discuss the azimuthal 
profile of the shock. 
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10%. The temperature of electrons is not as strongly af- 
fected by particle acceleration: with back-reaction, T c is 
lowered by about a factor of 2 in the shocked ISM (it re- 
mains the same within about 1% in the shocked ejecta). 
The different behaviour of electrons and protons with 
respect to particle acceleration stems from our hypothe- 
sis of progre ssive e quilibration downstream of the shock 
(see Section 2.1.2). Even though the ratio (3 = T e /T p 
ends up being much higher than the initial value m e /m p 
assumed at the shock, it is far from reaching unity, no 
matter how low T p gets: in the ISM it reaches about 1% 
in the un-modihed case and about 2% in the modified 
case. 

3.1.3. Ionization 

Knowledge of the density, density- weighted age, and 
electron temperature of the plasma allows us to compute 
the ionization sta te of the plasma a posteriori, as ex- 
plained in Section |2.1.3 On Figure [3] we show the total 
electron density. Comparing with the top of Figure [T 
it can be seen that in the ISM, which mostly consists of 
light, fully ionized, species, there is close to one free elec- 
tron per nucleon; whereas in the ejecta, there are many 
less free electrons, as heavy species take time to lose their 
electrons. We will not discuss in detail the precise mor- 
phology of all ionization states of all species, but instead 
present an average view of the state of a few important 
ones. On Figure H] we show maps of the mean electric 
charge < Z > of Oxygen, Silicon, and Iron. < Z > is de- 
fined as the sum of the charges of all the ionization states 
of a given element, weighted by their respective density 
fractions. It starts at < Z >= 1 at the shock fronts 
(by assumption), and quickly rises behind the shocks to 
reach a kind of plateau, at a value dependent on the total 
charge of the atom. The more electrons an element has, 
the less chances it has to approach full ionization. Super- 
imposed on this global radial evolution, Rayleigh- Taylor 
structures are clearly visible. 

The effect of acceleration on the mean electric charge 
mostly stems from the higher compression of the shocked 
ISM: with back-reaction, the region of highly ionized 
states starts closer to the forward shock, which is closer 
to the contact discontinuity. It is worth recalling that 
we only present here an average quantity < Z >, which 
can mask spatial variations in the fraction of different 
levels of ionization of a given element. Local variations 
in the plasma state, and the impact of acceleration, will 
become more apparent when looking at the thermal emis- 
sion pixel by pixel. 

3.2. Thermal X-ray emission 

Using the data shown in the previous section, the emis- 
sion code computes the thermal emission in each cell, 
which is then projected along the line of sight (the z- 
axis) to generate synthetic maps. This allows us, for the 
first time, to present realistic maps of the thermal emis- 
sion from a SNR, including the effects of hydrodynamic 
turbulence and of particle acceleration. We first show 
global maps an d spect ra over the broad energy band sim- 



ulated (Section 3.2. 1[ ); then we show the emission maps 
(and equivalent width maps when a pplicable) for several 
narrow energy bands (Section 3.2.2), chosen to illustrate 



the behaviour of some prominent emission lines and of 



the high-energy continua. Such maps can be compared 
with observations made with modern and future X-ray 
instruments allowing for spatially resolved spectroscopy. 

3.2.1. Broad band emission 

Maps. The broad-band emission from 300 eV to 
10 keV is shown on Figure p] (the emission is computed 
in the code from 100 eV, but below 300 eV interstel- 
lar absorption will prevent instruments to observe the 
theoretical fluxes and features simulated). The emission 
includes the emission of all species (from Hydrogen to 
Nickel) and in all media (both the ejecta and the ISM). 
At the top of the figure, we show the emission in slices 
(as in the previous sections), to show where it actually 
comes from inside the shocked region: mostly from the 
Rayleigh- Taylor fingers of ejecta, plus a small and rather 
uniform contribution of the ISM. At the bottom of the 
figure, we show projected maps (the novelty of this sec- 
tion) obtained by summing all the slices along the line 
of sight. A distinctive 'fleecy' structure is obtained in 
projection within the remnant, because of the Rayleigh- 
Taylor instabilities. Again, this demonstrates the need 
for 3D simulations to obtain realistic maps. We note that 
the outer edge of the ejecta is always brighter than the 
interior of the remnant, because more of the Rayleigh- 
Taylor structures is seen there in projection. After sum- 
mation over all the pixels, the total differential emissiv- 
ity of the octant of SNR simulated is of the order of 
10 39 ph.s _1 .eV _1 (that is a differential flux of the order 
of 10 -5 ph.cm _2 .s _1 .eV _1 at a distance of 1 kpc). Af- 
ter summation over the broad energy band considered, 
the total integral emissivity of the octant is of the order 
of 10 43 ph.s -1 (that is an integral flux of the order of 
10 _1 ph.cm _2 .s _1 at a distance of 1 kpc). 

The effect of acceleration can be seen on both the in- 
tensity levels and the patterns of the emissivity. The 
emissivity of the shocked ISM is always low (compared 
to the ejecta), with back-reaction at the forward shock 
it is even lower, plus the emitting region is substantially 
narrower. The emissivity of the ejecta remains about 
the same, as a result the contrast between the ejecta and 
the ISM becomes more pronounced around the contact 
discontinuity. And as the total emission is dominated 
by the ejecta, it is not much affected by particle back- 
reaction (by about 18% in the case simulated here). We 
also note that the general aspect of the ejecta is not visi- 
bly modified, although the Rayleigh- Taylor structures at 
small scales are all spatially re-distributed (a quantifica- 
tion of the degree of turbulence would be interesting, but 
is beyond the scope of this paper). 

Spectra by medium. To discuss more precisely the 
impact of acceleration on the ISM and the ejecta, on 
Figure [6] we plot the broad-band energy spectra spatially 
integrated over different regions. The top black curves 
show the emissivity of the entire region simulated - this 
offers an orthogonal view to Figure [5j instead of sum- 
ming all energy bins in each pixel, here we sum all pix- 
els in each energy bin. The blue (red) curves show the 
emissivity of a region dominated by the ISM (ejecta), de- 
fined as all the (projected) cells where the fraction of ISM 
(ejecta) is > 95%. We already noted that the emissiv- 
ity of the ISM is substantially lower than the emissivity 
of the ejecta. Now it becomes apparent that the spread 
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increases with energy. On the plot, we have indicated 
a few energy bands that are dominated by a particular 
emission line of a unique element - many more lines can 
be seen in the spectrum, but they blend with each other 
at the spectral resolution chosen. The emissivity of a 
line depends on the abundances of the elements, which 
are different in the ISM and in the ejecta (but which 
get mixed in the vicinity of the contact discontinuity). 
For instance, one can see prominent lines of Silicon and 
Sulfur in the ejecta. For a given element, the emissivity 
of the lines also depends on its ionization state, which in 
turn depends on the density, temperature, and age of the 
cell, and therefore varies strongly with position (see next 
section). As a result, spectra with distinctively different 
shapes (not shown here) can be obtained when extracting 
small spatial regions in the projected maps. 

The effect of acceleration is to reduce the intensity of 
the emission, all the more so at the highest photon en- 
ergies. The situation is very different for the ISM and 
for the ejectaj^] since we consider efficient acceleration 
to occur at the forward shock. For the ISM-dominated 
region, with back-reaction the emissivity is reduced by a 
factor of about 2 around 1 keV and by more than 2 or- 
ders of magnitude at 10 keV. This will have a strong 
impact on the visibility of the post-shock regions with 
an X-ray instrument. The total emissivity of the SNR is 
always dominated by the ejecta, though, and so is less 
affected by back-reaction: it is reduced by about 30% at 
100 eV and about 75% at 10 keV. As a result, the rel- 
ative visibility of the ISM and of the ejecta is strongly 
affected by particle acceleration: back-reaction enhances 
the contrast between the ejecta and the ISM. The dras- 
tic reduction of the emission of the ISM at high energies 
is explained by the exponential decrease of the contin- 
uum emission with temperature: the continuum emis- 
sion dominates at high energies, and the te mperat ure is 
much lower with back-reaction (see Section 3.1.2). We 



note that line emission is not affected in the same way: in 
fact, the few lines present at high energies become more 
prominent above the continuum. 

Spectra by Elements. To conclude this section on 
the broad-band emission, we decompose the total spec- 
trum of Figure [7] in another way, by separating the emis- 
sion of the various elements handled by the code. We 
do not show all 15 here, but only 3 of particular im- 
portance: Oxygen, Silicon, and Iron on Figure [7J As 
the atomic number increases, lines appear at higher and 
higher energies (provided the conditions are met in the 
plasma) . 

The effect of acceleration is fairly subtle when isolat- 
ing a single element in this way. The emission is always 
lower with back-reaction, but the differences are small, 
and mostly seen at the highest energies (at 10 keV, they 
reach 40-50%). We are showing here three heavy ele- 
ments, that are more abundant in the ejecta. Similar 
behaviour is obtained for the other elements. We must 
emphasize here that the spectra plotted on Figure [7] are 
spatially-integrated spectra, that show the global prop- 
erties of the remnant. Actually the emission varies a lot 

14 Here we warn the reader, before making direct comparisons 
between the spectra of the two media and in the two acceleration 
cases, that the number of extracted cells varies, as indicated in the 
legend of the figure. 



spatially because of the different conditions in density 
and temperature, and the back-reaction of particles does 
affect the local conditions, so that spectra for any ele- 
ment can have distinctively different shapes when look- 
ing at specific regions of the SNR. 

3.2.2. Emission in narrow energy bands 

As the former analysis showed that most features are 
energy-dependent, we now present emission maps for 
three of the narrow energy bands displayed on Figure [6] 
and [71 the 'O-K band' (545 - 580 eV) on Figure K| the 
'Si-tCband' (1715 - 1840 eV) on Figurej9) and the 'Fe- 



K band' (6250 - 6520 eV) on Figure 10 These bands 



were chosen from the presence of a prominent line of an 
element (for the conditions of this particular remnant at 
this particular age) and labeled with the name of this 
element, however they are defined by a range of energies, 
and will invariably contain a continuum emission from 
all elements. Because they are narrow bands centered at 
different energies, they probe different conditions of the 
plasma, and therefore different regions of the SNR, and 
different levels of the impact of particle acceleration. To 
make this even more apparent, we propose another kind 
of visualization on Figure 12 RGB composites made out 
of the emissivities in the U-K, Si-K and Fe-K bands (cor- 
responding to low, medium, and high photon energies, 
respectively). 

Line emissivity maps. At the top of the figure for 
each band, as well as on the composite figure, we show 
the total emissivity in the band (as on Figure [5]). Note 
that the intensity scale is different on each figure, be- 
cause it varies by 6 orders of magnitude. The turbulent 
pattern of the ejecta is visible on all emissivity maps, 
but with notable variations. The O-K and Si-K bands 
probe the deepest parts of the ejecta (i.e., the feet of 
Rayleigh- Taylor fingers) , while the Fe-K band probes the 
shallowest parts of the ejecta (i.e., the mushroom heads 
of Rayleigh- Taylor fingers). There is no radial strati- 
fication of the different elements in our initial setup, so 
what is observed here is the effect of the variations in the 
plasma conditions. This is caused by the fact that the Fe- 
K line requires high temperatures, which are found only 
close to the contact discontinuity, while the O-K and Si- 
K lines shine at pretty much all temperatures found in 
the shocked region, including close to the reverse shock 
(model emissivities in the four bands chosen, as a func- 
tion of temperature and ionization age, are provided for 
reference in Appendix [B]) . Hence observing the remnant 
at different energies, corresponding to different elements, 
allows to perform a radiography of its structure. 

The effect of acceleration is not to be seen on the global 
pattern of emission, but rather on the local variations in 
emission: the distribution of Rayleigh- Taylor structures 
is slightly modified with back-reaction. As before, the 
impact is more pronounced at high energies: variations 
are more visible in the Fe-K band. Because the tempera- 
ture is lower, with back-reaction the top of the Rayleigh- 
Taylor fingers is no longer bright, and the Fe-K map gets 
closer to the maps of other elements. These effects are 
not however very prominent on the emissivity maps be- 
cause the emission is proportional to the square of the 
density, and the density of the ejecta (which make most 
of the emission in terms of intensity and patterns) is not 
strongly affected by acceleration at the forward shock. 
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Equivalent width maps. To focus on the line emis- 
sivity, we now present maps of 'equivalent width' at the 
bottom of Figures J§J [9j [T0| These maps are produced by 
first subtracting the continuum emission from the total 
emission (in the energy band selected) to get only the 
emission of the line, then by dividing by the continuum 
emission to remove the strong dependence of the em ission 
on density ( [Hwang et al]|2000[ [Miceli et al.| |2006). The 
value obtained represents the width (in energy) of the un- 
derlying continuum that one needs to consider to obtain 
the same number of photons as in the line. The depen- 
dence on density being removed, this value is now only 
a function of the electronic temperature, of the abun- 
dance of elements (which we chose to be uniform in each 
medium), and on the ionization state of each element 
(which depends on density and temperature). It is re- 
markable that the maps obtained are not uniform: the 
ejecta pattern is still clearly visible (although it looks 
different than before). 

The effect of acceleration is much more visible on the 
equivalent width maps than on the emissivity maps. 
As before, back-reaction effects are more pronounced at 
higher energies: they are limited in the O-K band, but 
important in the Fe-K band. The effects are two-fold. 
First, the maximum width gets higher (by about 75% in 
the Fe-K band), which means that back-reaction tends 
to enhance the line emission with respect to the con- 
tinuum emission. This is caused by the fact that the 
line and continuum emission have a different dependence 
to the changes in temperature induced by particle ac- 
celeration. The main effect of back-reaction is to lower 
the temperature in the ISM, and therefore to lower its 
emission. Given that the ISM emission is mostly a con- 
tinuum emission, while most of the line emission comes 
from the ejecta, the result is that the diffuse ISM emis- 
sion is removed and the emission lines of the heavy ejecta 
become more evident. Secondly, the contrast decreases 
in the interior region of the SNR (by about 25% in the 
Fe-K band), which means that back-reaction tends to 
smoothen the line emissivity spatially. This is partly 
due to the fact that the shocked region shrinks, so that 
all features get compressed. 

Continuum emissivity maps. To conclude this sec- 
tion, we show on Figure emission maps in two broader 
energy bands, located athigh energies: the 'contl' band 
from 4 keV to 6 keV and the 'cont2' band from 8 keV 
to 10 keV. In our model these bands are devoid of line 
emission (note that we do not include the faint contribu- 
tions from Chrome and Manganese), and therefore illus- 
trate the fate of the continuum emission. The continuum 
emission directly scales with density, so that these maps 
are mostly maps of the (projected) density squared. The 
continuum emission also decreases exponentially at high 
temperatures, so that the maps are not exactly the same 
in the two energy bands. 

The effect of acceleration on the shocked ISM is high- 
lighted, because we have excluded emission lines which 
mostly trace the ejecta. As particles are accelerated at 
the forward shock, the temperature of the ISM is signif- 
icantly reduced downstream of shock, and thus its high 
energy emission. This can be seen on the projected maps 
when looking towards the edges of the remnant, between 
the contact discontinuity and the forward shock: without 
back-reaction the emissivity is low immediately down- 



stream of the shock and rises when getting closer to the 
ejecta, with back-reaction the shock front is closer from 
the ejecta, but the emissivity remains low everywhere 
downstream of it. This can also seen when looking to- 
wards the interior of the remnant (where the ISM is on 
top of the ejecta): with back-reaction, a veil of diffuse 
emission from the ISM is removed. 

4. DISCUSSION 

In this section we summarize the most salient features 
we have observed on the various simulated maps and 
spectra, and how they impact the interpretation of ob- 
servations. We also review the assumptions underlying 
our model, to assess how the results might differ with a 
different setup. 

4.1. General morphology and energy transfers 

Even before conducting a detailed computation of the 
thermal emission, one expects morphological modifica- 
tions to be visible, that are caused by the presence (and 
escape) of a significant population of relativistic parti- 
cles. In the simulation presented here, we estimate that 
particle acceleration took away about 20% of the kinetic 
energy released by the explosion. As the resulting fluid 
gets more compressibl e, the shocked region shrinks (De- 
courchelle et ah] 2000): it gets narrower and denser, as 
seen on Figure [I] The first effect: a more compact struc- 
ture can provide an observat ional diagnostic, as already 
shown in |Ferrand et al. ( 2010[ ) and associated references. 
The second effect: a higher compression ratio is more 
difficult to measure, because it requires the detection of 
the thermal emission from the plasma located just be- 
hind the shock, a region which was found to be usually 
dominated by non-thermal (synchrotron) radiation. 

Besides the density, the level of thermal emission 
critically depends on the temperature of the plasma, 
more precisely on the electrons temperature (T e ). The 
unavoidable effect of efficient particle acceleration at the 
shock is a decrease of the temperature of the shocked 
material. The drop is large for ions (Figure [2), which 
take part in the dynamics of the shock: in the case 
presented T p is reduced by up to one order of magnitude 
in the shocked ISM. However, the situation is different 
for electrons: there is no reason why they would be 
heated in the same way as protons at the shock, and 
their temperature will also depend on subsequent energy 
exchanges downstream of the shock. Assuming, as we 
do, a progressive Coulombian equilibration between 
electrons and protons, T e appears to be always limited 
by the age of the remnant. T e is still affected by particle 
acceleration, indirectly, through the decrease of T p : it 
is lowered by a factor of up to about 2 in the shocked 
ISM (Figure [2j) . If we were considering the limit case 
of instant equilibration of electrons with ions at the 
forward shock, then T e would be higher than in our 
model, but the effect of back-reaction would be even 
more drastic. We note that the obtained values of T e 
are in a ran ge w here emission lines can be very sensitive 
(see Figure |15[ ) , and that the continuum emission has 
an exponential dependence on T e , so that even small 
variations of T e can have a major impact on the thermal 
emission. 
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4.2. Energy- dependent effects: the need for spectroscopy 

We always observed that, as part of the energy is chan- 
neled into energetic particles at the shock, the thermal 
emission of the plasma is reduced. However, the overall 
impact is small (less than 10% in our simulation), and 
variations in the absolute value of the X-ray flux are un- 
likely to be used as signatures of particle acceleration 
given the uncertainties on all the relevant parameters. 
Much more reliable will be the comparison of fluxes in 
different energy bands. Indeed, we observed a strong en- 
ergy dependence of back-reaction effects, as can be seen 
on spatially-averaged spectra (Figure[6]) or by comparing 



maps in different energy bands (Figures |8j [9 10 11 1. As 
a rule, the higher the energy of the photons observed, 
the higher the reduction in the emissivity. As a result, 
observing the remnant at different energies probes differ- 
ent levels of impact of the back-reaction of particles on 
the shocked material. In general effects can be measured 
only above ~1 keV, and they are much more obvious for 
the ISM. Despite the increase in density downstream of 
the forward shock, we find that the high-energy emis- 
sion from this region is strongly reduced. This is due 
to the fact that the ISM emission is mostly a continuum 
emission, which has a strong dependence on temperature. 
Emission lines have a different dependence on tempera- 
ture, and so are not affected in the same way: actually 
the reduction of the continuum caused by back-reaction 
makes the lines more prominent (this can be seen on the 
spectrum of Figure [6] for the ISM as well as on the equiv- 
alent width maps of Figure 10 for the ejecta). 



4.3. Space- dependent effects: the need for 
spatially-resolved spectroscopy 

For a given element, the effects of back-reaction are 
barely visible on overall spectra (Figure [7| , where the 
emission is spatially-averaged over the whole SNR. But 
on the emission maps, spectra can look noticeably differ- 
ent from one pixel to another, and in the same location 
when considering back-reaction or not. In particular, it 
is important to separate the ISM and the ejecta, which 
have different emission properties, and are affected in a 
different way by back-reaction, but which are mixed up 
by Rayleigh- Taylor instabilities. Regarding simulations, 
this implies a 3D treatment to get the turbulent distri- 
bution of the shocked material. Regarding observations, 
this is a call for spatially-resolved spectroscopy. Given 
the strong variations in density and temperature that we 
observe on the maps, the average quantities that are com- 
monl}P°| derived from fits of the global spectra may be 
an oversimplification. For instance, in which sense the 
temperature derived from the spatially-averaged emis- 
sion can be regarded as a spatial average of the temper- 
ature of all the emitting cells should be clarified. Even 
when a 2D spatially-resolved study is conducted, we re- 
call the obvious fact the emission is integrated along the 
line-of-sight. Only 3D simulations can properly explain 
how the resulting features, as seen in projection with 
modern spectro- imaging instruments, are related to the 
local plasma conditions - as affected by particle acceler- 
ation. 

1 5 The observer can a lso perform multi-component fitting, see 
e.g. |Kosenko et al.| | (2uTo] . 



The spatial effects of back-reaction are much better 
seen on equivalent width maps rather than on total emis- 
sivity maps, after pure density effects have been removed, 
in order to reveal temperature and abundances effects 
on ionization and emission. We note that our synthetic 
equivalent width maps are not uniform at all, despite 
the use of uniform abundances (respectively in the ISM 
and in the ejecta): they map variations in the plasma 
conditions (temperature and ionization age), which have 
a non-trivial pattern in 3D (and even more when seen 
in projection) because of the hydrodynamic instabilities. 
So observing these could not be directly used to map 
the abundances of the elements, at least on small scales. 
If the various elements were not distributed evenly spa- 
tially, one could see bigger clumps of emission on the 
maps, as observed in some remnants. We deliberately 
kept the description of the ejecta simple in this paper, 
and regard the simulations presented here as a reference 
case. The next step will be to initialize the simulations 
with a non-uniform distribution of the ejecta, as obtained 
from supernova models, to reflect the combustion layers 
formed in the progeni tor stars and/or during th e explo- 
sion (see for instance Decourchelle et al. (2001) for ob- 
servational evidence for incomplete mixing ot the Silicon 
and Iron layers). This can be done in our current nu- 
merical setup by the addition of a passive tracer for each 
species we want to follow. 

4.4. The fate of the ejecta 

The impact of back-reaction of accelerated particles 
is distinctly different for the (shocked) ISM and for the 
(shocked) ejecta. This of course primarily stems from 
the fact that we consider acceleration at the forward 
shock only: the temperature of the ISM is much lowered, 
whereas the ejecta are only affected indirectly, through 
the dynamics around the contact discontinuity. Another 
aspect is that the two media are not affected in the same 
way by changes in the plasma conditions, because of their 
different composition: most emission lines of heavy ele- 
ments, present in the ejecta, are less affected by temper- 
ature than the continuum. As a result, the visibility of 
the ejecta with respect to the ISM is enhanced by the 
back-reaction. We note that indeed, on X-ray images 
of Tycho's SNR, the shocked ejecta shaped by Rayleigh- 
Taylor structures are clearly visible, whereas a shocked 
ISM c omponent is not req uired to explain the observa- 
tions ( | Warren et al. 2005 1 . On the other hand, on our 
simulated maps the contrast inside the ejecta seems to 
be somewhat reduced by the back-reaction (this is visible 
on the equivalent width maps of Figures P 



9 10) 



As the region between the contact discontinuity and 
the forward shock shrinks, the density and temperature 
profiles are modified on one side of the contact disconti- 
nuity, and the development of the Rayleig h- Taylor struc- 
tures can be affected. In our earlier work (Ferrand et al. 
2010 1 we did not observe a major impact on the over- 
all growth of the fingers. In this new study, looking at 
the X-ray emission in different energy bands allows to 
better see the effects of back-reaction on their structure. 
Radial variations of the line emissivities, caused by the 
evolution of temperature and ionization age downstream 
of the shocks, are visible on the RGB representations of 
Figure 12 With back-reaction, the low-energy emission 
in the O-K and Si-K bands gets more prominent and the 
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high-energy emission in the Fe-K band gets suppressed. 
So the morphology of the fingers is energy-dependent, 
and as a result is acceleration-dependent. We note that 
the Fe- K emission observed with Chandra in Tycho was 
used by |Warren et al. (2005) to trace the position of the 
reverse shock, wher eas m our synthetic maps the Fe-K 
emission (Figure 10 1 rather traces the location of the con- 
tact discontinuity, because only there sufficient electronic 



temperatures ar e reached for this high-energy line. War 
ren et al. ( 2005[ ) assumed the presence of hot Iron close 



to the reverse shock, based on the observation made with 
X MM- Newton, that the Fe-K emission peaks inside the 
Fe-L emission ( |Decourchelle et al.|[200T| ) . In our simula- 
tion, all Fe-L lines are faint. This is m agreement with the 
dependence on plasma conditions plotted on Figure |15| 
Fe-L line emission requires higher ionization ages that 
are difficult to reach in the ejecta. As a result, the Fe-L 
lines are always blended with stronger K lines from other 
elements (mostly Oxygen and Neon), except in a small 
energy banti around 730 — 750 eV (see spectra of Fig- 
ures m 



and [7]); and even there, the emissivity of the line 
is of"the same order as the emissivity of the continuum 
(mostly Hydrogen and Helium). Maps in this narrow 
energy band (not shown here) reveal that the Fe-L emis- 
sion can indeed extend further than the Fe-K emission, 
although at flux levels that are not observationally rel- 
evant. To have a prominent Fe-K line emission closer 
to the reverse shock, we would need to somehow reverse 
the temperature profile between the reverse shock and 
the contact discontinuity. This could be caused either 
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re reverse shock: the 



model of Tycho of Badenes et al. (|2005a|) works best for 
#3 - 0.03. 

Here we recall that a simplification made in our nu- 
merical model is to decouple the computation of the ion- 
ization of the plasma from the computation of its hydro- 
dynamic evolution: we record the evolution of the ion- 
ization age, and compute the actual ionization state at 
the end of the hydro simulation. The most accurate (but 
more numerically demanding) way to handle ionization 
would be to do the time integration of the balance equa- 
tions for each species concurrently with the time inte- 
gration of the conservation equations for the global fluid. 
A comparison of the two approaches can be found in |Pat-| 
naude et al. (|2010[ ) for the shocked ISM, based on ID 
simulations of a remnant in the Sedov phase. It is shown 
that post-processing the ionization leads to slightly over- 
estimating the final thermal X-ray emission, especially 
at low photon energies and for high plasma densities. 
We observed that back-reaction effects are the highest 
at the highest energies, and we chose a fairly small am- 
bient density, so we do not expect dramatic changes to 
our results regarding the relative effects of particle ac- 
celeration. We also note that the intensity of the ion- 
ization coupling depends on the fraction of massive ele- 
ments, so it is not expected to be critical for the shocked 
ISM (where back-reaction effects are the strongest in our 
model, by construction). In the ejecta the situation is 
more complicated, this is why we loop o yer th e compu- 
tation of the ionization (see end of Section|2.1.3 ) , in order 



particles at the reverse shock, the direct coupling of the 
ionization of at least the most abundant species should 
be investigated. 

4.5. Acceleration at the reverse shock? 

Acceleration of particles at the main blast wave of 
a SNR is well established, so it appears mandatory to 
explore the impact of the acceleration on the develop- 
ment of the shocked ISM. The situation is less clear for 
the reverse shock, that propagates backwards inside the 
ejecta. The major argument against considering accel- 
eration at the reverse shock is that, during the early 
expansion phase of the SNR, the magnetic field frozen 
in the ejecta gets diluted orders of magnitude below the 
levels required to accelerate particles. Magnetic turbu- 
lence is a key component of the diffusive shock acceler- 
ation mechanism: it traps particles close to the shock 
front, where they can gain energy. However, accelerated 
particles are believed to be able to generate themselves 
the magnetic turbulence that scatter them, by triggering 
various instabilities upstream of the shock. It is there- 
fore conceivable, if the acceleration process can start at 
a very early stage, that the magnetic field gets so much 
amplified that t he process is then se lf-sustained. We re- 
fer the reader to Ellison et al. (2005) for an investigation 
of efficient non-linear acceleration at the reverse shock, 



and to Telezhinsky et ah] (2012) for a recent study in 
the test-particle regime ( using ID hydrodynamic simu- 



lations for both works). |Telezhinsky et al" p012[ ) find 
that the reverse shock can make a significant contribu- 
tion of the total spectrum in the early expansion phase 
(up to about 400-600 years, which corresponds to our 
final age). Even in the linear regime, total spectra never 
look like single-shock test-particle solutions, so that the 
question arises of making the distinction between modifi- 
cations caused by non-linear effects at the forward shock 
and/or by an additional contribution from the reverse 
shock. Authors agree that the maximum energy reached 
by particles at the reverse shock is always lower than at 
the forward shock, so that we anticipate a lesser impact 
on the high-energy emission from the remnant. Still, re- 
cent radio and X-ray observations have been interpreted 
as the non-thermal emissi on from the reverse shock re- 
gion in som e SNRs: CasA dGotthelf et al. |2001 



Holder & 



Vink|2008|), Kepler peLahey ct al. 2002), RCWSBlRho 



to get consistent estimates of the final electron density 
and temperature. If we are to consider acceleration of 



et aj.|2 002)', IE 0102.2-7219 ( jSasaki et al.j2006| . We are 

not aware of any such indications for Tycho, which was 
used as our reference object. In any case, it seems ap- 
propriate to add acceleration at the reverse shock as a 
second step in future work. 

4.6. SN Type 

Finally, we recall that we considered in this paper the 
remnant of a fiducial type la supernova, hence the as- 
sumption of a uniform circumstellar medium. The as- 
sumption of perfect spherical symmetry might be the 
reason why we observe such a strong brightening at the 
edges of the shocked ejecta seen in projection, which is 
sharper than what is observed in remnants of the kind 
simulated such as Tycho's SNR. In reality, the medium in 
which the remnant expands will always exhibit some in- 
homogeneities. Such features will be specific to a particu- 
lar remnant though, whereas we are trying here to exhibit 
generic features of the acceleration and back-reaction. 
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In the case of the remnant of a core-collapse supernova, 
we would need to consider the radial density profile of the 
circumstellar medium, shaped by the wind of the mas- 
sive progenitor. This can be accounted for in our model, 
using a power-law profile of index s = 2 for the ISM. 
This substantially affects the hydrodynamic profiles in 
the shocked region, and therefore the thermal emission. 

5. CONCLUSION 

We have developed a 3D model that simulates the mor- 
phological and spectral evolution of a young SNR under- 
going efficient particle acceleration. To the best of our 
knowledge, this is the first time the thermal X-ray emis- 
sion is computed when considering conjointly (i) the 3D 
development of hydrodynamical instabilities (Rayleigh- 
Taylor) in the shocked region, (ii) the time-dependent 
back-reaction of particles accelerated at the blast wave, 
and (iii) the progressive temperature equilibration and 
non-equilibrium ionization of the plasma. 

Because we assumed efficient acceleration to take place 
at the forward shock only, back-reaction effects mostly 
impact the shocked ISM (an additional contribution on 
the shocked ejecta from the reverse shock might be useful 
to consider for some remnants). Our simulations repro- 
duce the most well-known effects: the shell of shocked 
material gets narrower, denser, and of lower tempera- 
ture. Note that the energetics of the blast wave directly 
affects the temperature of protons, whereas the emis- 
sivity of the hot plasma depends on the temperature of 
electrons. We assumed that electrons are not efficiently 
heated at the shock (which is supported by both obser- 
vations and other simulations, and can be regarded as 
a conservative hypothesis), and that electrons progres- 
sively equilibrate with protons through Coulombian in- 
teractions downstream of the shock (equilibrium is never 
reached at the age considered) . This complicates the dis- 
cussion of the impact of shock acceleration on the ther- 
mal emission. We observed that the electron tempera- 
ture in the shocked ISM is significantly reduced (here by 
a factor of about two). X-ray emission is still granted, 
but will be significantly affected by the efficiency of par- 
ticle acceleration. 

Because less energy goes into the thermal fluid, the ef- 
fect of particle acceleration, if any, is generally to lower 
the overall thermal emission. However, absolute varia- 
tions of fluxes are small, and unlikely to be used as diag- 
nostics given all the uncertainties in the modelling of any 
given object. More useful are the facts that variations in 
flux depend very much on the energy band considered, 
and that spectral features vary spatially. This demon- 
strates the need for spatially-resolved high-resolution 
spectroscopy, and the usefulness of such simulations. In 
our model of a young SNR evolving in a homogeneous 
ISM, we observed that, the higher the photon energy, the 
more pronounced the effects of back-reaction: they are 
barely visible below 1 keV, and are easily measurable at 
10 keV. The continuum and line fluxes are not affected 
in the same way, because of their different dependence 
on temperature. As a result, we observed that back- 
reaction enhances the visibility of the lines over the con- 
tinuum. Equivalent width maps (obtained by subtract- 
ing and dividing the continuum emission) were shown 
to be the most useful visualization of the emission, be- 
cause they probe the temperature and ionization age, and 



thus the level of back-reaction. Note that, even assuming 
uniform abundances, we obtain non-uniform equivalent 
width maps, a property that must be kept in mind when 
interpreting X-ray observations. 

On all sorts of maps, the interior region of the SNR 
bears the signature of the instabilities developing be- 
tween the shocked ISM and the shocked ejecta. Assum- 
ing the ejecta are of radially uniform abundances, the 
Ray leigh- Taylor fingers have a variable radial extent in 
X-rays, because of the global increase of temperature and 
ionization age towards the contact discontinuity, that re- 
flect the history of shock conditions and the transport 
downstream of the shock. As a result, their global mor- 
phology in narrow energy bands is affected by particle 
acceleration. A more qualitative study would be in order 
to assess how the small scale distribution of turbulent 
features is modified. Unsurprisingly, the level emissivity 
of the ejecta is not very much affected by efficient ac- 
celeration at the forward shock. However, the relative 
visibility of the ejecta is clearly enhanced by the back- 
reaction of particles on the shocked ISM, as the veil of 
diffuse emission from the ISM is removed. On the other 
hand, the contrast within the ejecta is reduced: back- 
reaction seems to smoothen interior regions of the maps. 

In this study, we used physical parameters most rel- 
evant for a type la SNR. However our model is fairly 
general, and could be adapted for a type II SNR. For 
simplicity we assumed uniform ejecta, and a uniform ISM 
(with different elemental abundances in each). This work 
should therefore serve as a reference case for future stud- 
ies where we will assess the impact (i) of a more complex 
distribution of the ejecta stemming from the explosion 
and (ii) of an imprint of the progenitor on the ambi- 
ent medium. Regarding the first aspect, |Orlando et al.| 



(20121 have recently shown that ejecta clumping has a 
stronger impact on the shock structure than efficient ac- 
celeration in SN 1006. Some of our present choices were 
inspired by Tycho's SNR, but we did not try to make 
an actual fit of the emission of any given SNR at this 
point. Actually, we do not reproduce the peak of Iron 
line emission close to the reverse shock, as suggested by 
X-ray observations of Tycho. This prompts a more de- 
tailed study of the evolution of the electron temperature 
in the shocked ejecta. The purpose of this paper was 
first to outline the method and demonstrate its useful- 
ness. We can now consider making a more precise model- 
ing of some well observed historical remnants. As such, 
the work presented here provides a framework for the 
interpretation of current SNR observations (with XMM- 
Newton, Chandra, Suzaku), as well as for the preparation 
of observations with future X-ray missions (in particular 
Astro-H). 

In this paper, we focused on the thermal emission from 
the SNR, as a way to diagnose the presence of energetic 
particles through back-reaction on the flow. The accel- 
eration model we are using also computes the spectra 
of particles accelerated at the shock. In a follow-up pa- 
per, we will present results on the non-thermal emission 
produced by these energetic particles present inside the 
shell. We will thus obtain a complete, broad-band, view 
of the morphology and X-ray emission from young SNRs. 
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APPENDIX 
A. EFFECT OF RESOLUTION 

We show here some maps with a coarser spatial grid, to investigate the role of numerical resolution. We first 
note that, regardless of the resolution, we can always make comparisons between the ejecta and the ISM, between 
cases with and without back-reaction of particles, and see the effects of progressive temperature equilibration and 
ionization downstream of the shocks. Numerical resolution is a concern mostly because of the fractal-like structure of 
the Rayleigh- Taylor instabilities: more and more details appear in the fingers as the grid cells become smaller. On 
Figure [l3j we show the same maps of the broad-band thermal emission as on Figure |5j with a grid of 256 3 instead 
of 1024 3 (that is 4x times less cell in each of the 3 spatial directions). Although the same overall structure can 
immediately be recognized, the two sets of maps look different: we can see finer structures at higher spatial resolution. 
However, the average spectra (not shown here) stay the same within a few %. Here a word of caution on the scales 
is in order: the maximum flux appears to be higher by more than one order of magnitude on Figure |l3"| compared to 
Figure [5] simply because of the lower pixel density (divided by 4 x 4 = 16). If we re-bin the high-resoTution maps by 

3lO( 



4x4 blocks to get the same pixels as on the low-resolution maps, then fluxes agree as can be seen on Figure 14 (fluxes 
in the re-binned slices are still 4x lower, because the slices are still 4x thinner). Remarkably, the smoothed maps 
from the high-resolution simulation contain much more details than the maps from the low-resolution simulation, even 
though both have the same effective pixel resolution. There are small-scale structures present in the high-resolution 
simulation, that do not appear in the low-resolution simulation even though there is room for them. This demonstrates 
that the growth of the Rayleigh- Taylor structures has entered a non-linear regime, where the development of structures 
at a given scale depends on the existence of structures at other scales. As a consequence, numerical simulations made 
to match the angular resolution of a particular instrument must actually be substantially over-sampled, in order to 
catch all the physics possibly observed. 

B. LINE EMISSIVITIES 

We show here for reference some line emissivities, computed from the thermal emission code of Section |2.1.4| as a 
function of temperature and ionization age. The emissivity is integrated over 4 different energy bands, corresponding 
to lines of Oxygen, Silicon and Iron. We recall that the electronic temperature T c in the shocked region for the remnant 
considered is of the order of 10 7 K: it raises from a few million degrees at the shock fonts to about 20 million degrees 
in the center. The ionization age goes up to n e t ~ 2.5 x 10 9 cm _3 .s. 

It is apparent that the emissivity depends very strongly on both parameters T e and n e t: it varies by ten orders 
of magnitude over the range considered. We observe the same distinctive curved shape for all elements, but with 
important differences. For instance, for Iron, triggering the Fe-K line requires high enough temperatures (T c > 10 7 K), 
but then is relatively weakly dependent of the ionization age (as long as n e t > 10 9 cm _3 .s); triggering the Fe-L line 
requires much stricter conditions (in particular T c must be very close to 10 7 K if n e t is in the range 10 9 — 10 10 cm~ 3 .s). 
Emission in the O-K and Si-K bands requires less restrictive conditions (it is granted over most of the range of T e and 
n e t we are concerned with). These behaviours help the interpretation of the RGB maps and equivalent width maps 
presented in Section [3.2. 2 1 
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Fig. 1. — Slices of hydro-dynamic quantities, in the 2 = plane, in the shocked region. Top: mass density, bottom: ionization age 
(time since the material was shocked, weighted by the electron density). Cases without ('OFF', on the left) and with ('ON', on the right) 
back-reaction of particles are compared. 




Fig. 2. — Slices of the temperature, in the 2 = plane, in the shocked region. Top: protons, bottom: electrons. Cases without ('OFF', 
on the left) and with ('ON', on the right) back-reaction of particles are compared. 




Fig. 3. — Slices of the electronic density, in the 2 = plane, in the shocked region. Cases without ('OFF', on the left) and with ('ON', 
on the right) back-reaction of particles are compared. 




Fig. 4. — Slices of the mean electric charge of a given element, in the 2 = plane, in the shocked region. Top: Iron, middle: Silicon, 
bottom: Oxygen. Cases without ('OFF', on the left) and with ('ON', on the right) back-reaction of particles are compared. 
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Fig. 5. — Total thermal emissivity in a broad 0.3 — 10 keV band. Top: slices in the 2 = plane, bottom: projected maps along the 
2— axis. On both types of maps, the differential emissivity at the source etot is plotted for each simulated pixel. The flux at the Earth will 
be diluted according to the square of the distance d to the source. The actual fluxes recorded by a given instrument will depend on its own 
angular resolution (we recall that the maps presented here have 1024 pixels along each direction, at the age shown here they span 5.61 pc, 
hence a physical cell size of ~ 0.0055 pc, and an angular resolution of about 1" at d = 1 kpc). Cases without ('OFF', on the left) and with 
('ON', on the right) back-reaction of particles are compared. 
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Fig. 6. — Broad-band energy spectra, spatially integrated over different regions, with or without the back-reaction of particles (solid or 
dashed lines). Top curves (in black) show the sum of the emissions from all cells (note that there are 12% less emitting cells in the ON 
case with respect to the OFF case). Middle curves (in red) show the contribution of cells dominated by the shocked ejecta, which make 
34% of all the emitting cells in the un-modified case (OFF), and 29% of the cells in the modified case (ON). Bottom curves (in blue) show 
the contribution of cells dominated by the shocked ISM, which make 13% of all the emitting cells in the un-modified case (OFF), and only 
6% of the cells in the modified case (ON). Some prominent lines are labeled by the emitting element for reference. 
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Fig. 7. — Broad-band energy spectra for three elements (in colours: Oxygen in red, Silicon in Green, Iron in Blue) and sum of the 
contributions of all elements considered (in black), with or without the back-reaction of particles (solid or dashed lines), integrated over 
the whole octant of the remnant simulated. Some prominent lines are labeled by the emitting element for reference. 




Fig. 8. — Projected maps, along the z— axis, of the thermal emission in the 'O-K band' (545 — 580 eV). Top: emissivity, bottom: ratio of 
lines over continuum emissivity (equivalent width maps). Cases without ('OFF', on the left) and with ('ON', on the right) back-reaction 
of particles are compared. 




Fig. 9. — Projected maps, along the 2— axis, of the thermal emission in the 'Si-K band' (1715 — 1840 eV). Top: emissivity, bottom: ratio 
of lines over continuum emissivity (equivalent width maps). Cases without ('OFF', on the left) and with ('ON', on the right) back-reaction 
of particles are compared. 




Fig. 10. — Projected maps, along the z— axis, of the thermal emission in the 'Fe-K band' (6250 — 6520 eV). Top: emissivity, bottom: ratio 
of lines over continuum emissivity (equivalent width maps). Cases without ('OFF', on the left) and with ('ON', on the right) back-reaction 
of particles are compared. 




Fig. 11. — Projected maps, along the z— axis, of the thermal emission flux in two high-energy continuum bands. Top: 8 — 10 keV, bottom: 
4 — 6 keV. Cases without ('OFF', on the left) and with ('ON', on the right) back-reaction of particles are compared. 
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Fig. 12. — RGB rendering of three line emissions. For each pixel, the value of the red / green / blue channel is assigned from the emissivity 
of resp. the O-K band / Si-K band / Fe-K band, as displayed on Figures fs| / [91 / 1 10| (linearly normalized to 256 levels). Regions of pure 
blue, for instance, are dominated by Fe-K line emission. Regions of yellow = red + green, are made out of a blend of O-K and Si-K line 
emission. Top: slices in the 2 = plane, bottom: projected maps along the z— axis. Cases without ('OFF', on the left) and with ('ON', on 
the right) back-reaction of particles are compared. 




Fig. 13. — Total thermal emissivity. Same as Figure[5] from a simulation with spatial resolution 4 times lower in each direction. 




Fig. 14. — Total thermal emissivity. Same as Figure [3] re-binned by 4x4 blocks. 
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Fig. 15. — Line emissivities in four energy bands: O-K band: 545 - 580 eV, Fe-L band: 730 - 750 eV, Si-K band: 1715 - 1840 eV, 
Fe-K band: 6250 — 6520 eV. The quantity plotted still has to be multiplied by the density squared to get the actual emissivity per unit 
volume, hence the units of photons. cm 3 . s — 1 . 



